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We show that the increments of generalized Wiener process, useful to describe non- 
Gaussian white noise sources, have the properties of infinitely divisible random processes. 
Using functional approach and the new correlation formula for non-Gaussian white noise 
we derive directly from Langcvin equation, with such a random source, the Kolmogorov's 
equation for Markovian non-Gaussian process. From this equation we obtain the Fokker- 
Planck equation for nonlinear system driven by white Gaussian noise, the Kolmogorov- 
Feller equation for discontinuous Markovian processes, and the fractional Fokker-Planck 
equation for anomalous diffusion. The stationary probability distributions for some 
simple cases of anomalous diffusion are derived. 
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1. Introduction 

Stochastic dynamics is useful to model many biological, chemical, economical and 
physical systems. The random driving forces have very different origins, in most 
cases they are Gaussian white or colored noise sources, but often these forces must be 
considered as non-Gaussian ones, for example, in sensory and biological systems [1]. 
Moreover, in many physical and biological systems a deviation of real statistics of 
fluctuations from Gaussian law, leading to anomalous diffusion, is observed [2,3]. A 
suitable mathematical model to describe the anomalous diffusion is the fractional 
Fokker-Planck equation, which can be derived from different theoretical approaches 
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[3-7]. 



In this paper we show that the increments of generalized Wiener process, rep- 
resenting the integral of non-Gaussian white noise, have the properties of infinitely 
divisible random processes. This enables us to obtain the characteristic functional 
of non-Gaussian white noise. Then, by using functional approach to split the corre- 
lation between stochastic functionals, we derive the Kolmogorov's equation directly 
from Langevin equation with non-Gaussian white noise source. From this general 
equation we can obtain the Fokker-Planck equation for continuous Markovian pro- 
cesses, the Kolmogorov-Feller equation for discontinuous Markovian processes and 
the fractional Fokker-Planck equation for anomalous diffusion. 

2. Generalized Wiener process 

It is well known that Wiener process rj (t) has stationary and independent incre- 
ments, and its time derivative £ (t) = fj (t) is Gaussian white noise. The correlation 
function of the increments rj (t + T) — rj (i) depends on time difference T only, and 
increments of non-overlapping time intervals are statistically independent. Let us 
consider now a non-Gaussian random process with stationary and independent in- 
crements which can be called as generalized Wiener process. We can represent the 
increment t]t = rj (T) — 77 (0) as 



where At = T/n. Since the random variable t\t can be divided into the sum 
of arbitrary number of independent and identically distributed random variables, 
its probability distribution belongs to the class of infinitely divisible distributions. 
Hence, we can express the logarithm of characteristic function of the random value 
r/T in the Levy-Khinchine form [8] 



where Mx{dx} is the canonical measure which, according to definition, takes posi- 
tive values on a finite intervals, and the integrals 



converge for any positive x and T. Introducing the non-negative density Pt{x) of 
the canonical measure: Mx{dx} = px(x)dx, defined in generalized sense, we can 
rewrite Eq. J2J) as 



n 




(1) 







For two consecutive time intervals T and S, we have 



9 T +s (u) = (e 



< 



e m(jj T +jjs) 



} = (e iu ^)(e iur > s ) = e T (u)6 8 (u) , 
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and 

<t>T+S (u) = <\>t (u) + <ps (u) • (4) 
Because u is an arbitrary parameter, we obtain from Eqs. © and 

Pt+s (x) = pt (x) + ps (x) , Vx. (5) 

The differentiable solution of Eq. JSJ, regarding parameter T, is only linear one 

Pt (x) = Tp (x) . 

So, from Eq. we obtain 

<p T (u)=T - 2 p{x)dx, (6) 



where the kernel /? (x) > 0. 

3. Functional correlation formula for non-Gaussian white noise 

The time derivative of generalized Wiener process £ (t) = fj (t) is a stationary ran- 
dom process and has analogy to Gaussian white noise. Now we derive the char- 
acteristic functional of this non-Gaussian delta-correlated noise. By definition, we 
have 



Ot [u] = yexp |i J £ (t) u (t) dr jy = yexp |i J u(r)d77(r) 

= /expjt^lm ^(^[^(TibJ-f/CTTb-!)]!^ (7) 

= , Um n ( IT ex P { iu (^ fc ) ^ ( T k)-V ( T fc-i)]} ) = lim TT °Ar k (u (i?fc)) , 
\fc=l / fe=l 

where is some internal point of time interval (rfc_i,Tfe) , <5 T = maxAr^, Ar^ = 

Tfe ^ 7fc-i (to = 0, r„ = t), and ^Ar fc is the characteristic function of increments. 
To obtain Eq. (JZJ) we used the statistical independence of increments of generalized 
Wiener process ij (t). Further from Eqs. © and 10 we obtain 

^ri ,. tt f* r°° e iu(ti k )x _ i _ j u sin a; 
t [u\ = hm ||exp^Ar fe / - p(x)dx 

{« roc e tu(# h )x _ I _ fa g J nx "I 
lim > Ar fe / p (x) dx > 
6 ^°fi J-°° x J 

f /•* Z" 00 e te W x -l-m(r)sinx .. ,1 
= cxpl dr j — ^ p (x) dx ^ . (8) 



Now we derive a useful correlation formula for stochastic functionals of non-Gaussian 
white noise. We use the generalization of Furutsu-Novikov formula [9] , to split the 
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correlation between a stochastic functional R t [£] of non-Gaussian random process 
£ (r), denned on the observation interval r 6 (0,t), and the process £ (i) itself 



(£(*)i? t K + z]> = 



*t [u] 



■'(*) 



(fit K + *]>. 



(9) 



Here z (t) is an arbitrary deterministic function, and $t [u] = In 0j [u]. From Eq. (jSJ) 
we obtain the following expression for variational operator in Eq. JHJ 



e iu(t)a; _ j _ j M sin , 
(t) x 2 



■p (x) dx 



p{x) 



dx 



Substituting this expression in Eq. @ we arrive at 



p(x) 



dx 



exp< y 



Sz (t) 



cos y 



[ e *«(*)w -cosy]dy. 



(R t [£ + z])dy. 



By inserting the operator of functional differentiation into the average and by 
putting z = 0, we get finally 



it®* [£]} 



cxp < y 



^ K] ) - (Rt K])cosy 



(10) 



4. Derivation of Kolmogorov's equation 

Let us consider now the Langevin equation with a non-Gaussian white noise £(t) 

x = /(x,t)+.g( 2 ;,i)£(i). (11) 

By differentiating with respect to time the following expression for probability den- 
sity of the random process x(t) 

W(x,t) = (6(x-x(t))), (12) 

and taking into account Eq. (|11|1 . we obtain 

dW d d 

= \f(x,t)W] -q^9 (x, t) (£ (t) S(x-x (t))> . (13) 

By using functional differentiation rules and following the same procedure of Ref. [10] , 
from Eq. (|ll|l we get 



S d 
-^Jjj S(x-x (t)) = -— g (x, t) S (x - x (*)) . 



(14) 



Thus, the operator d/8£(t) of functional differentiation with respect to the func- 
tion 5 (x — x (i)) is equivalent to the ordinary differential operator —d/dx [g (x,t) ). 
Using Eq. H10JI in Eq. 1)13(1 and taking into account Eq. I|14|l . we obtain the Kol- 
mogorov's equation for nonlinear system l|ll|l driven by non-Gaussian white noise 



dW _ d[f(x,t)W] 
dt dx 



d 



d 



exp <j -z-g^g {x,t) \ -1 + sin z—g (x, t) 



dx- 



dzW, 
(15) 
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which is the main result of our paper. By series expansion of the exponential 
operator in Eq. I|15(l . we have 



dW 
~dt 



n—1 



d_ 

dx 



g(x,t) 



W. 



where 



At 



z — sm z 



p(z) dz, 



z n ~ 2 p(z)dz, n>2. 



(16) 



(17) 



The equation Qlfijl is valid if all integrals in Eq. Ijl7(l converge. It coincides with 
the result obtained in Ref. [10] for Poissonian white noise £(t), if we interpret the 
coefficients A n as the moments of random jumps of x(t). Moreover, this equation 
can be recast in the Kramers-Moyal form 



~dt ~ ^ 



(-if d n 

n\ dx r ' 



[K n (x, t) W] , 



(18) 



where K n (x,t), n = 1,2,... arc the kinetic coefficients. 



5. Particular cases of non-Gaussian noise sources 

Following Ref. [8] we analyze different kernel functions p(x) to obtain particular 
cases of Kolmogorov's equation l(T5)l. related to different non-Gaussian noise sources. 

(a) As a first simple case we consider a Gaussian white noise £(t). The cor- 
responding kernel function is p (x) — 2D5 (x) . Therefore we obtain the ordinary 
Fokker-Planck equation 



dW d d d 

— = -- l , {x ,t ) W ] + D- g{ x,t } - g ( 1 :,t ) W. 

(b) For Poisson kernel p (x) = B 8 (x — q), Eq. (|15f) becomes 



(19) 



dW 
~dt 



d_ 

dx 



f(x,t) + 



Bg (x,t) (q - sing) 



(-l) n q ri 



n=1 



dx 



g(x,t) 



W. 



(20) 



(c) When the increment r)x of generalized Wiener process has the gamma dis- 
tribution with exponent p, the kernel is p(x) — pxe~ ax (x > 0). As a result, after 
integration we find from Eq. (|15fl 



dW d 

— = -— [/ (x, t) - pg (x, t) arc cot a]W - p In 



18.' 
a ox 



W. (21) 



(d) For one-side a-stable Levy distributions the kernel function is p (x) = pg x 1 a 
(x > 0), and we get 

dW d 

[f{x,t)W] 



+Po 



dz 



at 

d 



dx 



d 



exp \ -z— g (x,t) } - 1 + — g (x, t) sin z 



dx 



W(x,t). 
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After rearrangements we obtain 

dW d[f(x,t)W} , po d f°° dy 

g{x,t) 



dt dx a dx ' J y 



w. 



(22) 

This equation is the fractional Fokker- Planck equation, and describes anomalous 
diffusion in the form of asymmetric Levy flights. 

(e) For additive driving noise £ (t), g(x,t) = 1 in Eq. and the exponential 
operator in Eq. Ijl5(l reduces to the space shift operator. As a result, we find 



W(x-z,t)-W (x, t) + sinz dW (X ' *' 



dz. 



dx 

"(23) 

Equation 1231) is similar to the Kolmogorov- Feller equation for purely discontinuous 
Markovian processes [4, 11] 

dW f°° 

— =v Q{x-z)W{z,t)dz-vW{x,t), (24) 

J — oo 

where Q (x) is probability density of jumps step, and v is the constant mean rate 
of jumps. Putting / (x,t) — and comparing Eq. ill'ol) with Eq. I|24|) we obtain the 
kernel function for this case 

p(x) = vx 2 Q{x). (25) 

For non-Gaussian driving force ^ (t) , with symmetric a-stable Levy distribution, 
the kernel function is p(x) = -D|ir| 1-Q . As a result, Eq. ll'MI) takes the following 
form 

and describes the anomalous diffusion in form of symmetric Levy flights. 

(f ) Finally, for multiplicative noise £ (t) with g (x, t) — x, the operator in Eq. i|15|) 
reduces to a scaling operator, and the generalized Kolmogorov's equation becomes 

d ~aT = ~Yx [f {X ' t] ~ KX] W + f T ^ W {e ~ ZXl £) ~ W {X -' t)] dZ ' (2?) 
where 

K = r p{z) t nz dz. 



6. Simple examples of superdiffusion 

Now we derive the stationary probability distribution Woo (x) for some simple cases 
of anomalous diffusion in a fixed potential U(x) and in overdamped regime (see 
also [6, 12, 13]). The Langevin equation for the particle displacement x(t) reads 

x = -U'(x)+t(t), (28) 

where £ (t) is a symmetric a-stable Levy process. The corresponding Kolmogorov's 
equation is Eq. H2(i(l . with f(x,t) — —U'(x). The stationary probability distribution 
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Woo (x) can be found from the following equation 

d Itt> ( \vir i . n r W °° ^ - W °° ( x ) a n roo^ 
— [U (x) W^+D ■ dz = 0. (29) 

J-oo |X — -2| 

a) For uniform potential well with two reflecting boundaries at x — ±L we get, 
from Eq. (|29|l . the uniform stationary probability density 

W^x) = ± 

as for ordinary Brownian diffusion. 

b) For linear system with parabolic potential U (x) = fix 2 /2, by Fourier trans- 
form of Eq. (|29|) . we obtain the asymptotic characteristic function of the process 
x(t) as [6] 

^W=exp{- 2r '-°)j'W2) |tr |. (30) 

c) For symmetric potential U(x) = j\x\, after Fourier transform of Eq. (|29|l we 
find 

Trey 

jkG (k) + 2D \k\ a T (-a) F (k) cos — = 0, (31) 

where 

/•oo />oo 

F(k)= / VFoo (a;) cos fcx dx, G(k)= I W 00 {x)sm.kxdx. 
Jo Jo 

Functions F (k) and G (k) are connected by Hilbert formula 

G(k) = -v.p. (°° f^dg, (32) 
7T 7 ft - g 

where w.p. means the principal value of integral. After substitution of Eq. I|32|) in 
Eq. I|31|) we obtain finally the following integral equation for cosine Fourier transform 
of stationary probability density 

v-p. r P^- dq + ^-\k\ a T (-a) F (k) cos ^ = . (33) 
J k - q -fk 2 

d) For quartic potential U (x) — /3x 4 /4, from Eq. (|29l) we obtain the following 
differential equation for stationary characteristic function 

(3k^^-2D\k\ a T(-a)9 oo (k)cos^ = 0. (34) 
Equations (|33|l and (|34(1 can be solved numerically. 
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7. Conclusions 

Starting from the generalized Wiener process, whose increments have the probability 
distribution belonging to the class of infinitely divisible distributions, we obtained 
the new correlation formula for stochastic functionals. This allows us to derive 
the Kolmogorov's equation (|15(l for diffusion, induced by non-Gaussian white noise 
source, directly from Langevin equation. The Fokker-Planck equation for ordinary 
Brownian motion, the Kolmogorov-Feller equation for purely discontinuous Marko- 
vian processes and the fractional Fokker-Planck equation for anomalous diffusion 
are obtained as particular cases of the Kolmogorov's equation (|15|l . This equation 
can be used for determining the stationary probability distributions, probability 
distributions of first-passage times and other statistical characteristics of nonlinear 
dynamical systems. 
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